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ABSTRACT 

Aims. The aim of this paper is to study with high angular resolution a dense core associated with a low-luminosity IRAS source, 
IRAS 00213+6530, in order to investigate whether low mass star formation is really taking place in isolation. 

Methods. We carried out observations at 1.2 mm with the IRAM30m telescope, and VLA observations in the continuum mode at 
6 cm, 3.6 cm, 1.3 cm and 7 mm, together with HtO maser and NH3 lines toward IRAS 00213+6530. Additionally, we observed the 
CCS Jf^ = 2i-lo transition, and H2O maser emission using the NASA 70 m antenna. We studied the nature of the centimeter and 
millimeter emission of the young stellar objects (YSOs) found in the region, and the physical properties of the dense gas and dust 
emission. 

Results. The centimeter and millimeter continuum emission, together with the near infrared data from the 2MASS allowed us to 
identify three YSOs, IRS 1, VLA 8A, and VLA 8B, with different radio and infrared properties, and which seem to be in different 
evolutionary stages. IRS 1, detected only in the infrared, is in the more advanced stage. On the other hand, VLA 8A, bright at centime- 
ter and millimeter wavelengths, coincides with a near infrared 2MASS source, whereas VLA 8B has no infrared emission associated 
with it and is in the earliest evolutionary stage. The overall structure of the NH3 emission consists of three clouds. Two of these, MMl 
and MM2, are associated with dust emission at millimeter wavelengths, while the southern cloud is only detected in NH3. The YSOs 
are embedded in MMl, where we found evidence of line broadening and temperature enhancements. On the other hand, the southern 
cloud and MM2 appear to be quiescent and starless. Concerning the 1.2 mm dust emission, we modeled the radial intensity profile of 
MMl. The model fits reasonably well the data, but it underestimates the intensity at small projected distances from the 1.2 mm peak, 
probably due to the presence of multiple YSOs embedded in the dusty envelope. There is a strong differentiation in the relative NH3 
abundance with low values of ~ 2 x 10"** toward MMl, which harbors the YSOs, and high values, up to 10"', toward the southern 
cloud and MM2, suggesting that these clouds could be in fact in a young evolutionary stage. 

Conclusions. IRAS 00213+6530 is harboring a multiple system of low-mass protostars, indicating that star formation in this cloud is 
taking place in groups or clusters, rather than in isolation. The low-mass YSOs found in IRAS 00213+6530 are in different evolution- 
ary stages suggesting that star formation is taking place in different episodes. 
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1. Introduction ters (e. g.. lciarke et alj2000HLada & Ladal2003h including low 



mass stars (e. g.,|G6 mez et al. 1993; Huard et al.| ll999l : lLee et alj 

It is widely accepted that there are two modes of star formation: 20061: iBrookeetd.. .2007: Teixeira et al. 2007.). ind icating that 

the isolated mode and the clustered mode. This classification truly isolated star formation is rare.|Adarns&Myers|(|Mi pro- 

results from studying the association betweei i dense cores and pose an intermediate case between the isolated and clustered 

young stellar objects (YSOs). For example, [ Benson & Myers| ^odes, i.e., star formation in groups, and propose that most stars 

09891) study a wide sample of ammonia cores and its relation gj-Qupg ^nd/or clusters, i.e., in cluster environments, 

with the position of IRAS sources, and find that typically only However, the theories of low mass star format ion assum e that 

one IRAS source is associated with a single ammonia core in star fo rmation takes place in the isolated mode (IShu et al.ll 19871 

the Taurus Molecular Cloud, where star formation can be as- Lid^^99). Since star formation in cluster environments may 

sumed to take place in isolation. On the contrary Orion and jiff^er from the isolated mode (e.g., Pfalzneret al. 2008), it is 

Perseus would be examples of molecular clouds forming stars necessary to study with high angular resolution dense cores as- 

in clustered mode (|Lada etalj[l993). A broad base of recent sociated with one single IRAS source to assess if star forma- 

studies, carried out with higher angular resolution than that of jjon really taking place in isolation. In this context, we aim 

Benson & Myers, show that most stars form in groups or clus- investigating with high angular resolution a dense ammo- 

nia core associated with a single low luminosity IRAS source, 

Send offprint requests to: Gemma Busquet, IRAS 00213+6530 (hereafter 100213). 
e-mail: gbusquet@am.ub.es 
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Table 1. VLA observational parameters in the IRAS 00213+6530 region 
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100213, with a luminosity of < 20 L0 and at 85 pc of dis- 
tance , belongs to the molecular cloud M120. 1+3.0 dYang et al.l 
I1990I) in the CepheusOB4 star-forming reg ion. The region i s 
physically related to the Hii region S171 dYanget al.lll990h . 
The NH3 emission in the north of M120. 1+3.0 was stud- 
ied through single-dish observations by lSepulvedal (l200lh . The 
NH3 emission consists of two condensations, each one peak- 
ing very close to the position of an IRAS source, 100213 and 
IRAS 00217+6533 (100217), suggesting that both IRAS sources 
are deeply embedded in high density gas. The mass derived for 
the condensation associated with 100213 is > 45 Mq. The am- 
monia emission engulfing both IRAS sources is associated with 
CO high-velocity emission , indicating the presence of a molec- 
ular outflow in the region dYang etalJll990l) . However, it is not 
clear which IRAS source is driving the outflow. 

In this paper we report on high angular resolution observa- 
tions with the Very Large Array (VLA) of the continuum emis- 
sion at 6 cm, 3.6 cm, 1.3 cm, and 7 mm, as well as of the dense 
gas traced by NH3 (1,1) and NH3 (2,2) together with observa- 
tions of H2O maser emission. In addition we also present the 
continuum emission at 1 .2 mm observed with the IRAM 30 m 
telescope, and CCS and H2O maser observations carried out with 
the NASA 70 m antenna at Robledo de Chavela. The paper lay- 
out is the following: in § 2 we describe our observations and the 
data reduction process, and present the main results for the con- 
tinuum and molecular line emission in §3. In §4 we analyze 
the dust and NH3 emission and show the method used to derive 
the NH3 abundance in this region. Finally, in § 5 we discuss our 
findings, and we list the main conclusions in § 6. 

2. Observations 

2.1. IRAM 30 m observations 

The MPIfR 37-element bolometer array MAMBO at the 
IRAM 30 m telescop43 was used to map the 1 .2 mm dust con- 
tinuum emission toward 100213. The observations were carried 
out in 2006 June 2. The main beam has a HPBW of 10". We 
used the on-the-fly mapping mode, in which the telescope scans 
continuously in azimuth along each row. The sampled area was 
200" X 140", and the scanning speed was 5" s"'. Each scan 
was separated by 5" in elevation. The secondary mirror was 
wobbling at a rate of 2 Hz in azimuth with a wobbler throw 
of 46". The average zenith opacity was in the range 0.3-0.4. 
Pointing and focus were done on NGC 7538. The rms of the fi- 

' IRAM is supported by INSU/CNRS (France), MPG (Germany), 
and IGN (Spain) 



nal map was ~ 3.8 mJy beam . Data reduction was performed 
with the MOPSICJI software package that contains the necessary 
scripts for data reduction (distributed by R. Zylcka). 

2.2. VLA radio continuum observations 

The observations were carried out using the VLA of the NRACEI 
in the D configuration in the continuum mode at 6 cm and 3.6 cm 
on 2000 September 23, and at 3.6 cm and 7 mm during 2004 
August 24. The observational parameters for each epoch are 
summarized in Table [1] During the first epoch absolute flux cal- 
ibration was achieved by observing 3C286, with an adopted flux 
density of 7.49 Jy at 6 cm and 5.18 Jy at 3.6 cm. The absolute 
flux calibrator was 00137+331 (3C48) during the 2004 obser- 
vations, for which flux densities of 3.15 Jy and 0.53 Jy were 
adopted at 3.6 cm and 7 mm, respectively. In order to minimize 
the effects of atmospheric p hase fluctuations, at 7 mm we used 
the fast switching technique (ICarilli & Holdawavll997.) between 
the source and the phase calibrator over a cycle of 80 seconds, 
with 50 seconds spent on the source and 30 seconds on the cal- 
ibrator The 1.3 cm continuum emission was observed on 2006 
December 2 simultaneously with the H2O maser emission (see 
below. Sect. 2.3). 

Calibration and data reduction were performed using stan- 
dard procedures of the Astronomical Imaging Processing 
System (AlPSflof the NRAO. Clean maps at 3.6 cm and 7 mm 
were mad e using t he ta sk IMAGR of AIPS with the robust pa- 
rameter of'Briggs' ('1995) set equal to 5, which is close to natural 
weighting, whereas the map at 6 cm was made with the robust 
parameter equal to zero in order to obtain a synthesized beam 
similar to that at 3.6 cm. Since the signal-to-noise ratio of the 
longest baselines of the 7 mm data was low, we applied a uv- 
taper function of 80 kA in order to improve the sensitivity. 

2.3. VLA NH3 and H2O maser observations 

The observations of {J,K) - (1,1) and (J,K) = (2,2) inver- 
sion lines of the ammonia molecule were carried out in the same 
run as the 2004 continuum observations. In Table [T] we sum- 
marize the observational parameters. The adopted flux density 
of the absolute flux calibrator 0137+331 (3C48) was 1.05 Jy at 

- See http://www.iram.es/IRAMES7mainWiE7CookbookMopsic 
' The Very Large Array (VLA) is operated by the National Radio 
Astronomy Observatory (NRAO), a facility of the National Science 
Foundation operated under cooperative agreement by Associated 
Universities, Inc. 
* See h ttp://aips.nrao.edul 
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Table 2. Parameters of the 1.2 mm emission 



Position 






a (J2000) 6 (J2000) 


J peak 




(h m s) (° ' ") 


(mjy beam"') 


(mJy) 


MMl 00 24 10.7 65 47 09 


146 


332±35 


MM2 00 24 04.4 65 48 12 


46 


146± 9 



a wavelength of 1.3 cm, and the bandpass calibrator used was 
0319+415 (3C84). We used the 4IF spectral Hne mode, which al- 
lows simultaneous observations of the NH3 (1,1) and (2, 2) lines 
with two polarizations for each line. The bandwidth used was 
3.12 MHz, with 63 channels with a channel spacing of 48.8 kHz 
(0.6 km s ' at 1.3 cm) centered at vlsr - -19.0 km s plus a 
continuum channel that contains the average of the central 75 % 
of the bandwidth. 

The water maser line at 22.2351 GHz (6i6-523 transition) 
was observed with the VLA in the C configuration during 2006 
December 2. The phase center was the same as for the NH3 ob- 
servations, and the adopted flux density of the absolute flux cali- 
brator, 0137+331 (3C48), was 1.13 Jy at 1.3 cm. We summarized 
the observational parameters in Table [T] We used the 4IF mode, 
employing two IF with a total bandwidth of 3.12 MHz, with 63 
channels with a channel spacing of 48.8 kHz (0.6 km s ') cen- 
tered at -10 km s"' 0, and two IF with a total bandwidth of 
25 MHz were used to observe the continuum emission. 

The NH3 and H2O data were reduced with the standard AlPS 
procedures. The images were constructed using natural weight- 
ing in both cases. 

2.4. NASA 70m CCS and H2O maser observations 

We carried out a H2O maser emission monitoring toward 
IRAS 00213+6530 with the NASA 70m antenna (DSS-63) at 
Robledo de Chavela (Spain). The observations were performed 
in 2008 Apr 18, June 19, and September 23 using a cooled high- 
electron-mobility transistor as 1.3 cm front-end, and a 384 chan- 
nel spectrometer as backend, covering a bandwidth of 16 MHz 
(~ 216 km s"' with 0.6 km s"' resolution). Spectra were taken 
in position-switching mode. The HPBW of the telescope at this 
frequency is ~ 41". The typical system temperature was 120 K 
and the total integration time was around 30 min (on+off) for 
each session. 

In addition to the H2O maser observations, we also observed 
the CCS Jn - 2i-lo ti-ansition (22.344 GHz), with the same 
bandwidth and spectral resolution used for H2O. The CCS tran- 
sition was observed on October 3rd 2008 and on February 7th 
2009 during a total integration time of 26 and 40 minutes, re- 
spectively. The system temperature was 60 K and 73 K, respec- 
tively. 

For all the observations, the rms pointing accuracy of the 
telescope was better than 10". A noise diode was used to cali- 
brate the data, and the uncertainty in the flux calibration is es- 
timated to be ~ 30 %. The data reduction was performed using 
the CLASS package, which is part of the G1LDAS0 software. 



^ Note that the center velocity for the H2O maser line observed with 
the VLA is shifted by 9 km s~' from th a t of th e NH3 observations since 
the HtO maser reported by iHan et 2] 0191) was detected at vlsr = 
-0.7 km S-' 

^ See ,http://www.iram.fr/IRAMFR/GILDAS] 
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0'^4"'20' 24'°15' 24°'10' 24"05" 24°'00" 
Right Ascension (JSOOO) 

Fig. 1. 1.2 mm continuum emission toward IRAS 00213+6530. 
Contour levels are 3, 6, 9, 12, 15, 20, 25, 30 and 35 times the 
rms of the map, 3.8 mJy beam The synthesized beam of the 
IRAM30m is shown in the bottom left corner of the image. 
White crosses indicate the position of the millimeter sources 
VLA 8A and VLA 8B, and the white tilted cross marks the posi- 
tion of IRS 1. VLA 7 is indicated by the black cross (see Sect. 3.2 
for the objects labeled in this figure). 



3. Results 

3. 1. Continuum at 1.2 mm 

Figure [T] shows the 1 .2 mm continuum emission observed with 
the IRAM30 m telescope toward 100213. The overall structure 
of the dust emission consists of a central and compact dust con- 
densation, MMl, with some extended structure to the west. A 
2D Gaussian fit to MMl yields a deconvolved size of 13" x 12", 
P.A.= 92°. In addition, we also detected a fainter dust condensa- 
tion, MM2, located to the northwest of MMl, and elongated in 
the north-south direction, in a filamentary structure connecting 
MMl and MM2. The main results are summarized in Tabled 

3.2. VLA radio continuum emission 

We detected 11 sources at 6 cm, and 7 sources at 3.6 cm above 
the 4 cr detection threshold. Figure|2]shows the 6 cm and 3.6 cm 
continuum emission observed with the VLA toward 100213. In 
Table|3]we list the positions and flux densities, corrected for pri- 
mary beam response, of the detected sources, and the estimated 
spectral indices between 6 cm and 3.6 cm. To obtain the 3.6 cm 
map we analyzed separately the observations of the two epochs 
(see Sect. 2.2) in order to see the degree of variability of the 
detected sources at this wavelength. The final map was obtained 
after subtracting the uv data of VLA 3, which presents a variabil- 
ity larger than 8cr, and VLA 6, detected in 2004 but not during 
the observations carried out in 2000. In Table|4]we show the flux 
density measured in 2000 and 2004 for VLA 3 and VLA 6, as 
well as the variability during this period. Both sources, VLA 3 
and VLA 6, have a negative spectral index (from simultaneous 
observations at 6 cm and 3.6 cm) and are probably non-thermal 
extragalactic background sources. 
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Table 3. Parameters of the continuum sources detected in the IRAS 00213+6530 region 





Position" 
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^ Positions taken from the 3.6 cm map, except for non-detected sources at this band, for which positions correspond to the 6 cm map. 
Source well outside the primary beam at this wavelength. Primary beam correction larger than 12, and the uncertainty in the corrected flux is 
very large. 

The upper limit for non-detected sources is 4cr. 

Flux density measured in 2000 September 23. Source highly variable, so to compute the spectral index we used observations carried out at the 
same epoch. 
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Fig. 2. In both panels grey contours represent the main beam brightness temperature of the main line of the NH3 {J,K) - (1,1) 
inversion transition from Sepiilveda (2001). Left: VLA 6 cm continuum emission map (black contours) of the 1002 13 region. Contour 
levels are -3, 3, 4, 6, 10, 14, 18, 22, 32, and 42 times the rms of the map, 24 /iJy beam"^ The synthesized beam, 17'.'9 x ir.'7, 
with P.A.= 54°, is shown in the bottom left corner of the image. Right: VLA 3.6 cm continuum emission map (black contours) of 
the 100213 region. Contour levels are -3, 3, 4, 6, 8, 10, and 12 times the rms of the map, 18 /iJy beam"'. The synthesized beam, 
15'.'14 X 9'.'89, with P.A.= -72°, is shown in the bottom left corner of the image. The positions of the IRAS sources are indicated 
by crosses. IRAS 00213+6530 is at the center of the figure, while IRAS 00217+6533 is at the north-east. The VLA primary beam, 
9' at 6 cm and 5' at 3.6 cm are also indicated by black circumferences. 



As can be seen in Fig. [3] (middle panel), at 3.6 cm we de- 
tected two sources toward 100213, VLA 7 and VLA 8, separated 
by 20" and both inside the position error ellipse of the IRAS 
source. These sources are barely resolved at 6 cm, with VLA 8 
just being a weak prolongation to the north-east of VLA 7 (see 
Fig.[3]top panel). Yet, the higher angular resolution of the 3.6 cm 



map allows to separate the two sources. VLA 8 peaks close to 
the position of the dust condensation MMl, whereas VLA 7 lies 
~ 20" to the south-west, in the extended structure of the dust 
emission (see Fig. [T]i- Both sources are spatially resolved and 
VLA 8 shows a weak tail extending to the west. A 2D Gaussian 
fit to the two sources (excluding the weak tail of VLA 8) yields 
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Table 4. Highly variable sources at 3.6 cm 
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Fig. 3. Top: VLA contour map of the 6 cm continuum emission 
toward 100213. Contour levels are -3, 3, 6, 9, 12, 15, 18, and 21 
times the rms of the map, 24 /iJy beam"'. Middle: VLA 3.6 cm 
continuum emission. Contour levels are -3, 3, 4, 5, 6, 7, 8 and 9 
times the rms of the map, 1 8 fiiy beam" ' . The synthesized beams 
at 6 cm (17'.'9x 1 1'.'7, with P.A.= 54°), and at 3.6 cm (15'.'1 x9'.'9, 
with P.A.= -72°), are shown at the bottom left corner of each 
panel. Bottom: White: VLA 7 mm continuum emission, black: 
VLA 1.3 cm continuum emission. Contour levels at 7 mm are 
-3, 3, 4, 5, and 6 times the rms of the map, 0.2 mJy beam"'. At 
1.3 cm contour levels are -3, 3, 4, 5, 6, and 7 times the rms of 
the map, 0. 15 mJy beam"' . The synthesized beams at 1 .3 cm and 
7 mm are shown at the bottom left and right corners of the im- 
age, respectively. In all panels grey scale is the 2MASS TiTs-band 
image, and the position error of IRAS 00213+6530 is indicated 
by the ellipse. 



deconvolved sizes of 9'.'7 x 2'.'2 (P.A.= 46°), and 13'.'4 x 3'.'1 
(P.A.= 63°) for VLA 7 and VLA 8, respectively. The spectral 
index in the 6-3.6 cm range of VLA 7 is -1.6±0.2, characteris- 
tic of non-thermal emission, whereas VLA 8 has a spectral in- 
dex of 0.5 ±0.4, which is consistent with free-free thermal emis- 
sion from ionized gas that may be arising from a thermal radio 
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Variability 
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(mJy) 
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VLA 3 


0.38 ± 0.07 


1.86 ±0.07 


79% (~ 2I0-) 


VLA 6 


< 0.25" 


0.41 ± 0.05 


>39%(> 3a-) 



The upper limit for non-detected sources is 4(X. 



jet. Positive spectral indices, i.e., a > -0.1, have been found 
to be associated with sources driving m olecular outflows (e. g., 
lAnglada et al.]ll998l: iBeltran et alJl2001h . The 2MASS /T^-band 
image shows two sources, 2MASS J00241 1 10+6547095 and 
2MASS J00241010+6547091, the first nearly coinciding with 
VLA 8 and the second, named IRS 1, lying ~ 6" to the west. 

The maps of the 1.3 cm and 7 mm continuum emission 
obtained with natural weighting are shown in Fig. [3] (bottom 
panel). While we did not detect 7 mm continuum emission to- 
ward VLA 7, the 7 mm emission of VLA 8 is resolved into 
two components, VLA 8A and VLA 8B, separated by ~ 5" 
(~ 4300 AU at the distance of the source). At 1.3 cm we de- 
tected one source associated with VLA 8B, whose peak posi- 
tion coincides within ~ 0'.'5 with the 7 mm peak (Fig. [3] bottom 
panel). The 1.3 cm source is elongated roughly in the northeast- 
southwest direction, and is spatially resolved only in one direc- 
tion, with a deconvolved size of 1'.'2 (~ 1020 AU), at PA.= 21°. 
In Table |5] we show the position, peak intensity, and flux density 
of VLA 8 A and VLA 8B at 1.3 cm and 7 mm, as well as the 
spectral index between these wavelengths. In order to properly 
estimate the spectral index in the 1 .3 cm-7 mm range we applied 
a uv-taper function of 80 kA to obtain similar angular resolu- 
tions at both wavelengths (see Sect. 2.2). The resulting spectral 
indices are > 1.4 and 0.2±0.6 for VLA 8 A and VLA 8B, respec- 
tively. 

As seen in Fig. [3] (bottom panel), VLA 8A coincides with 
flie near infrared source 2MASS J00241 1 10+6547095, whereas 
VLA 8B has no infrared emission associated with it. In addition, 
~ 6" (~5100 AU ) west of VLA 8A there is the near-infrared 
source IRS 1 with no detected 7 mm emission. At 7 mm both 
thermal dust emission and free-free emission can contribute to 
the total emission. Then, in order to estimate the mass of gas and 
dust from the 7 mm emission, we need to know first the flux den- 
sity coming from thermal dust emission. Since the angular reso- 
lution achieved at 3.6 cm is not sufficient to resolve the two mil- 
limeter sources we cannot estimate the contribution of free-free 
emission at 7 mm for each individual source. A first approach 
is to smooth the 7 mm emission to the angular resolution of the 
3.6 cm data. The result is a single source with a flux density of 
2.5 ± 0.9 mJy, which is consistent with the sum of the flux den- 
sities of both millimeter sources. Extrapolating the flux obtained 
at 3.6 cm for VLA 8 to millimeter wavelengths with the spectral 
index obtained from the centimeter emission at 6 cm and 3.6 cm 
(a=0.5), we find that the expected free-free emission at 7 mm is 
~ 0.48 mJy (19 % of the 7 mm flux), indicating that the thermal 
dust component dominates at this wavelength. We note that the 
free-free contribution estimated at 7 mm can be considered as an 
upper limit because the free-free spectral index between 3.6 cm 
and 7 mm is expected to be flatter than between 6 and 3.6 cm. 
However, the thermal dust contribution at 7 mm may be different 
for VLA 8A and VLA 8B. The fact that VLA 8A is not detected 
at 1.3 cm above a 4 cr level indicates that at 7 mm the emis- 
sion is mainly due to thermal dust emission rather than free-free 
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Table 5. Parameters of sources detected at 1.3 cm and 7 mm in the IRAS 00213+6530 region 





Position" 


1.3 cm 


7 mm 




Spectral 




a (J2000) 6 (J2000) 


5, 


J peaK 


5,, 


Index*^ 


Source 


(h m s) (° ' ") 


(mjybeam"') (mjy) 


(mJy beam"') 


(mJy) 


1.3 cm-7mm 


VLA8A 


00 24 11.01 65 47 09.9 


< 0.6" 


1.22±0.26 


1.57±0.54 


> 1.4 


VLA8B 


00 24 11.76 65 47 09.8 


1.09±0.13 1.25±0.25 


1.12±0.38 


1.43±0.79 


0.2±0.6 



" From the 7 mm map. 
"4 0^ upper limit. 

° Estimated from the 1.3 cm flux density obtained with a uv-taper function of 80 k/i (see main text) 
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Fig, 4, Top panel: VLA channel maps of the NH3 (1,1) main hne. Contours levels are -6, -3, 3, 6, 9, 12, 18, 24, 30, 33, and 36 
times the rms noise of the map, 1.1 mJy beam Bottom panel: VLA channel maps of the NH3 (2,2) main line. Contour levels are 
-3, 3, 6, 9, and 12 times the rms noise of the map, 1.1 mJy beam In both panels the synthesized beam is shown in the bottom left 
corner. The positions of VLA 7 and VLA 8 are indicated by crosses. 



emission, whereas VLA 8B, which is associated with a 1.3 cm 
source, must have less thermal dust emission associated with it. 

To estimate the mass, we assumed that the dust emis- 
sion is optically thin, and used the opacity law of / Cy - 
0.01(v/230 GHz/ cm^g-i (lOssenkopf & HenningI [T99l . ex- 
trapolated to 7 mm. We used a dust emissivity index p - \.5 
(derived from the spectral energy distribution, see Sect. 4.2). The 
dust temperature is estimated by correcting the rotational tem- 



perature derived from NH3 (~ 20 K, see next sectio ns) to kinetic 
temperature (~ 25 K), following the expression of ' Tafalla et al.l 
(2004). Using the fraction of the 7 mm flux density arising from 
thermal dust emission (~ 2.1 mJy), the total mass derived for the 
two sources, VLA 8A and VLA 8B, is 6.3 Mq. It is worth noting 
that this mass is an upper limit since at 7 mm we are sensitive to 
spatial scales of ~ 2500 AU, smaller than that achieved with the 
NH3 observations. Thus, the temperature should be higher than 



G. Busquet et al.: Star formation in IRAS 00213+6530 



7 




O 



I I I I I 




I I I I I 



2 K 



I I I I I I 



IRSl 



MM2 



I I I I I I 



MM1 




sc 



- IK 




-30 -20 



-10 -30 -20 -10 
Velocity (km/s) 



Fig. 5. Spectra toward six positions of the IRAS 002 13 -1-6530 re- 
gion for NH3 (1,1) {left) and NH3 (2,2) (right), averaged over 
one beam. The six positions, which are labeled on the right panel 
of each row, are, from top to bottom, VLA 8 A, VLA 8B, IRS 1, 
MM2 (peak of the northwestern cloud), MMl (peak of the cen- 
tral cloud), and SC (peak of the southern cloud). The vertical 
scale for each transition is indicated in the bottom row. 



that estimated from the NH3 emission, and the resulting mass 
would be lower. The uncertainty in the mass is around a factor 
of 4 mainly due to uncertainties in the dust opacity and the dust 
emissivity index. 

As interferometers are not sensitive to large-scale struc- 
tures, we compared the 7 mm continuum emission with the 
1.2 mm dust emission from t he IRAM30 m tele scope. We fol- 
lowed the method described in lGirart et al.l (l2000h . which relates 
the FWHM of the single-dish emission with the half -power (u, v) 
radius, to estimate t he magnitude of this effect at 7 mm. By ap- 
plying the relation of lOirart et alj JIOO O), and adopting a FWHM 
for MMl of 0FWHM ^ (13" x 10")'^^ =11" the half-power (u, v) 



65°4S'00" - 



65°47'00' 



65°46'00' 



65°48'00' 



65°47'00" - 
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Fig. 6. Top panel: Zero-order moment of the NH3 (1,1) main 
line emission. Contours start at 1 %, increasing in steps of 10 % 
of the peak intensity, 0.0508 Jy beam"' km s"'. Bottom panel: 
Zero-order moment map of the NH3 (2, 2) main line emission 
(black contours). Contours start at 1 %, increasing in steps of 
15 % of the peak intensity, 0.0151 Jy beam ' km s"'. Grey con- 
tours: 1.2 mm continuum emission. Contour levels are 3, 6, 9, 
12, 25, and 35 times the rms of the map, 3.8 mJy beam"'. In 
both panels white crosses mark the position of the two millimeter 
sources, VLA 8A and VLA 8B, and the white tilted cross marks 
the position of the infrared source IRS 1 . The synthesized beams 
for each transition are shown at the bottom right corner, and the 
synthesized beam of the 1 .2 mm continuum data is shown in the 
bottom left corner of the bottom panel. 



radius of the 1.2 mm dust emission MMl becomes r ^ 8 kA. 
Given the shortest baseline of the VLA in the D configuration, 
which is 2.5 k/l < 8 kA, and the size of the observed emission 
one can estimate the fraction of correlated flux detected by the 
interferometer. For a source of 11", this corresponds to 94 % for 
our VLA configuration at 7 mm, indicating that at 7 mm we are 
not filtering out too much dust emission. 



3.3. NH3(1,1) and (2,2) 

NH3 (1, 1) and NH3 (2, 2) emission is detected in individual ve- 
locity channels from -19.6 to -20.9 km s"' and from -20.2 to 
-20.9 km s"', respectively. We also detected the inner satellite 
lines as well as one of the outer satellite lines of the (1,1) transi- 
tion, and we marginally detected one of the inner satellite lines 
of the NH3 (2,2). Figure ID shows the velocity channel maps of 
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Table 6. NH3 line parameters obtained from the Gaussian fits to the NH3 (1,1) and (2,2) 



Position Main line Inner satellites 





a (J2000) 


6 (J2000) 


Transition 


V'LSR 




7'mb 


Av 


Tun 


Av 


'MB ("!) 
^MB ('.') 


Source 


(h m s) 


(° ' ") 




(km s" 


') 


(K) 


(km s-') 


(K) 


(km s"') 


VLA 8A 


00 24 11.01 


+65 47 09.9 


(1,1) 
(2,2) 


-19.98 + 
-19.92 + 


0.01 
0.03 


3.9+0.2 
1.8+0.1 


1.19 + 0.02 
0.83 + 0.07 


1.8+0.2 


1.14+0.04 


2.2 


VLA SB 


00 24 11.76 


+65 47 09.8 


(1,1) 
(2,2) 


-19.89 + 
-19.88 + 


0.01 
0.04 


3.6+0.2 
1.5+0.1 


1.16 + 0.02 
0.91 + 0.08 


1.5+0.2 


1.18+0.06 


2.4 


IRS 1 


00 24 10.10 


+65 47 09.1 


(1,1) 
(2,2) 


-19.97 + 
-19.83 + 


0.02 
0.05 


2.2+0.1 
1.0+0.1 


1.19 + 0.03 
0.92 + 0.08 


1.0+0.1 


1.11+0.07 


2.2 


MM2 


00 24 12.30 


+65 46 45.0 


(1,1) 
(2,2) 


-19.84 + 
-19.85 + 


0.02 
0.08 


1.4+0.1 
0.5+0.1 


1.19 + 0.05 
0.88 + 0.17 


0.9+0.1 


1.06+0.05 


1.6 


MMl 


00 24 10.80 


+65 47 06.0 


(1,1) 
(2,2) 


-19.94 + 
-19.96 + 


0.01 
0.04 


4.8+0.2 
1.8+0.1 


1.10 + 0.03 
0.83 + 0.10 


2.5+0.1 


1.02+0.04 


1.9 


SC 


00 24 05.30 


+65 47 47.0 


(1,1) 


-19.95 + 


0.01 


2.7+0.1 


0.97 + 0.02 


1.5+0.1 


0.95+0.05 


1.8 




00^24^22f 18' 16' 14' 12' 1 O' 08' 06' 04' 02' 
a(J2000) 

Fig. 7. Top panel: First-order moment map of the NH3 (1,1) 
main line emission. Bottom panel: Second-order moment map 
of the NH3 (1,1) main line emission. Symbols are the same as 
in Fig.|6] Color wedge scales are km s The synthesized beam 
is shown at the bottom right corner of the images. Note that the 
second-order moment gives the velocity dispersion, and must be 
multiplied by the factor 2 V21n2 ^ 2.35 to convert to full width 
to half maximum. 



the NH3 (1,1) and (2,2) main line emission, and in Fig. |5]we 
show the NH3 (1,1) and (2, 2) spectra, not corrected for the pri- 
mary beam response, at some selected positions. In Table |6] we 
show the line parameters toward these positions obtained from a 
Gaussian fit to the NH3 (1,1) main and inner satellites lines, as 



Table 7. NH3 line parameters from the fits to the NH3 (1,1) mag- 
netic hyperfine components 





V'LSR 




Av" 






Source 


(km s" 


') 


(km s-') 


(K) 




VLA 8 A 


-19.99 + 


0.01 


0.8 + 0.1 


11.04+0.32 


2.71+0.13 


VLA8B 


-19.92 + 


0.01 


0.9 + 0.1 


8.41+0.35 


2.12+0.16 


IRS 1 


-19.92 + 


0.01 


0.9 + 0.1 


8.39+0.34 


2.11+0.15 


MM2 


-19.85 + 


0.01 


0.7 + 0.1 


6.60+0.34 


4.40+0.37 


MMl 


-19.96 + 


0.01 


0.7 + 0.1 


16.34+0.54 


3.46+0.17 


SC 


-19.96 + 


0.01 


0.6 + 0.1 


10.08+0.49 


3.96+0.27 



" Intrinsic line width (FWHM) of the magnetic hyperfine component. 
A = f(Jy(T^^) - /y(rbg)), where / is the filling factor, is the 
excitation temperature of the transition, T],^ is the background radi- 
ation temperature, and Jy(T) is the intensity in units of temperature, 
UT) = (hv/k)/(exp(hv/kT) - 1). 
Optical depth of the main line obtained from the fit. 



well as the NH3 (2,2) line. We additionally show the ratio of the 
main line to the inner satellites, which gives an indication of the 
line optical depth. The values obtained for the optical depth are 
in the range t,„ ~ 2 - 3. In Table |7] we present the line parameters 
resulting from the hyperfine fit to the NH3 (1,1) line toward the 
same positions. Given our spectral resolution, Av =0.6 km s"', 
the intrinsic line width has been obtained by selecting the value 
that minimize the hyperfine fit residual. The optical depths de- 
rived from the magnetic hyperfine fit are systematically higher 
but compatible with the values obtained from the Gaussian fit 
(see Appendix for further details). 

In Fig. |6] we present the zero-order moment map of the 
NH3(1,1) and (2,2) emission (integrated intensity including 
only the main line). While the NH3(1,1) map shows signifi- 
cant extended structure, the emission of the NH3 (2, 2) is com- 
pact around the position of the centimeter/millimeter sources. 
We also detected faint NH3 (2,2) emission toward the northern 
peak of MM2. The overall structure of the NH3 (1, 1) emission 
consists of three clouds. The central NH3 cloud coincides with 
the dust condensation MMl. Around ~ 22" south of MMl there 
is another cloud detected only in NH3 (hereafter southern cloud). 
Moreover, the NH3 cloud located to the northwest of MMl is as- 
sociated with the dust condensation MM2. Note that, because of 
the small size of the VLA primary beam at 1.3 cm, the sensitiv- 
ity toward MM2 is two times lower than toward the center of the 
field. Then, the dense gas emission traced by the NH3 molecule 
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Reference 


Observation 
date 


Telescope 


Beam 


Positional 
accuracy 


Detection 


rms 
(Jy beam"') 


-^peak 

(Jy beam"') 


Felli et al. 1992 


1991 Jan 18 


32 m Medicina 


1:9 


15" 


no 


1.1 


< 5.5 


Anglada et al. 1997 


1990 Feb 


37 m Haystack 


I'A 


15" 


no 


1.6 


<4.8 


Han et al. 1998 


1993 Nov 21 


13.7 m Purple Mount 


4.2' 


20" 


yes" 


9.1 


38.8 


this work 


2006 Dec 2 


VLA 


2'.'1 


< 1" 


no 


0.003 


< 0.015 


this work 


2008 Apr 18 


NASA 70 m 


40" 


10" 


yes'' 


0.03 


0.111 


this work 


2008 Jun 19 


NASA 70 m 


40" 


10" 


no 


0.04 


< 0.21 


this work 


2008 Sep 23 


NASA 70 m 


40" 


10" 


no 


0.06 


< 0.28 



" Detection at -0.7 km s ' . 

'' Marginal detection at -15 km s 



0.10- HjO 6,5-5. 




-100 -so 

Velocity (km s ' ) 



Fig. 8. Spectrum of the H2O maser observed with the 
NASA 70 m telescope on 2008 April 18. The spectrum was 
smoothed to a spectral resolution of 1.13 km s"'. The dashed 
line is a Gaussian fit to the spectrum. 



roughly follows the 1.2 mm dust continuum emission, except 
in the southern cloud in which we did not detect dust emission 
at all (see Fig.|6]). While the NH3 cloud MMl engulfs the cen- 
timeter source VLA 8 (i.e., IRS 1, VLA 8A, and VLA 8B), 
no NH3 emission is seen toward the position of VLA 7, which 
falls close to the edge of the NH3 emission. It is important to 
note that there is a near-infrared source detected in the 2MASS 
bands, 2MASS J0024 125 1+65464 18, spatially coinciding with 
the southern NH3 cloud, which probably is not associated with 
the dense gas (see Sect. 5.2 for a complete explanation). 

In Fig. I2] (top) we show the first-order moment (intensity 
weighted mean vlsr) of the NH3 (1,1) main line emission. As 
can be seen in this figure, the velocity along the central am- 
monia cloud MMl shows only small variations, and no sig- 
nificant velocity gradients are found between MMl and the 
southern cloud. The millimeter source VLA 8B is redshifted by 
~ 0.3 km s~' with respect to VLA 8A. In addition toward the 
western and eastern edges of MMl the NH3 emission is red- 
shifted by ~ 0.4 km s"'. Toward MM2 there is a small velocity 
gradient in the north-south direction of ~ 0.6 km s"' along a 
region of ~ 30". 

The map of the second-order moment (intensity weighted ve- 
locity dispersion) of the NH3 (1,1) main line emission is shown 
in Fig. I2] (bottom). The typical value found for the velocity dis- 
persion in the southern cloud is ~ 0.25-0.3 km s"', which cor- 
responds for a Gaussian line profile to a full width at half max- 
imum (FWHM) of 0.6 km s"', similar to the instrumental reso- 
lution. In contrast, in the NH3 cloud MMl there is evidence of 
line broadening toward the three embedded sources (forming an 



arc-shaped structure), with values up to 0.4-0.5 km s , corre- 
sponding to line widths of 0.7-1 km s"', corrected for instru- 
mental resolution, significantly higher than the expected thermal 
line width ~ 0.23 km s"' (estimated for a kinetic temperature 
of ~ 20 K), indicative of a significant contribution from non- 
thermal processes. We found that the typical value of velocity 
dispersion in the NH3 cloud MM2 is ~ 0.3-0.35 km s"', cor- 
responding to line widths of 0.4-0.5 km s"\ corrected for in- 
strumental resolution, and rises up to 0.45 km s"' (line width 
of 0.9 km s"', corrected for instrumental resolution) toward the 
northern peak of MM2. These values are slightly higher than 
those found in the southern cloud, suggesting that the gas is be- 
ing perturbed in the NH3 cloud MM2. 

Finally, we compared the VLA NH3(1,1) emission with 
the single-dish NH3 (1,1) emission of Sepiilveda (2001). From 
the VLA NH3 observations we estimated a peak intensity of 
27.8 mJy beam"', which corresponds to a main beam brightness 
temperature Tmb - 4 K. The largest angular scale detectable by 
the VLA at 1.3 cm in the D configuration is around 60", and the 
size of the largest features detected by us with the VLA is ~ 30". 
From the size of the NH3 emission we can estimate the dilution 
effect when observed this emission with a single-dish telescope 
of 1 'A. The dilution effect would decrease the main beam bright- 
ness temperature, Tmb. by a fact or of 7.8, consistent with the 
Tmb measured bv lSepulvedal (12001) . indicating that the fraction 
of emission filtered out by the interferometer must be small or 
negligible. 

3.4. H2O maser emission 

In Table [8] we compiled the H2O maser observations carried out 
toward 100213 up to date, including our observations with the 
VLA and NASA 70 m. The H2O maser em ission was clearly de- 
tected only in 1993 bv iHan et alj (Il998l) . using the 13.7 m ra- 
dio telescope of Purple Mountain Observatory, with a peak in- 
tensity of 38.8 Jy. We did not detect H2O maser emission to- 
ward 100213, except in the NASA 70 m observations on 2008 
April 18, where we marginally detected emission at a 3 cr level 
of 0.1 Jy beam"' (see Fig. The integrated intensity was 
~ 0.49 Jy beam"' km s"', and the velocity of the feature was 
-15.2 km s"', offset by ~ 5 km s"' from the velocity of the 

cloud, 20 km s"'. 

The characterization of H2O maser emission is difficult be- 
cause of its high temporal variability, and thus it is possible that 
there was no maser during the epochs of observation (2006 Dec 
and 2008) with the exception of the marginal detection during 
2008 April. However, it is n ot clear whether the H2O maser 
detected bv lHan et al] (Il998l) is associated with the 100213 re- 
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Fig. 9. Spectrum of the CCS Jn - 2i-lo transition observed 
with the NASA 70m telescope toward 100213. The dashed line 
is a Gaussian fit to the spectrum. 



Table 9. Parameters of the NASA 70 m CCS line emission to- 
ward 100213 



Line 


vlsr 




Av 


Tub 




(km 


') 


(km s-') 


(K) 


CCS = 2]— Iq 


-20.6 ± 


0.2 


1.3 ±0.5 


0.22±0.05 



gion. On one hand, the maser emission was detected at vlsr= 
-0.7 km s"', 19.2 km s"' offset from the systemic velocity of 
the cloud studied in this work. On the other hand, the beam of 
the 13.7 m radio telescope at this wavelength is 4'2, with a point- 
ing accuracy of 20", making it difficult to ascertain whether the 
maser is associated with the IRAS source. In addition, other at- 
tempts to detect H 2O maser emission toward 100213 have failed 
jFelli et al.ll 19921; C. Codella observed the source between 1989 
and 1999 using the 32m Medicina telescope but did not detect 
it; private communication). Thus, the maser activity of 100213 
seems to be currently in a rather quiescent state. 

3.5. CCS emission 

In Table |9] we show the parameters obtained from a Gaussian fit 
of the CCS line detected with the NASA 70 m, after combining 
the data of the two days observed, and we show its spectrum in 
Fig- HI The line is centered around -20 km s the same veloc- 
ity as NH3, and its line width (FWHM) is ~ 1.3 km s larger 
than the expected thermal line width for a kinetic temperature 
of ~ 20 K (estimated from NH3), which is of ~ 0.13 km s"'. 
This could be indicative of the CCS line width having a strong 
contribution from non-thermal processes, such as turbulence in- 
jected by the passage of an outflow and/or global systematic mo- 
tions. The measured line width is higher than the largest line 
width measured by de Gregorio-Monsalvo et al. (2006) toward 
a sample of low-mass YSOs. A high angular resolution study 
of Ide Gregorio-Monsalvo et al.l (12005) in CCS shows that this 
molecule is possibly enhanced via shocked-induced chemistry, 
and has a velocity gradient in the same direction of the outflow. 

4. Analysis 

4.1. Rotational temperature and column density maps 

We computed maps of the rotational temperature and column 
density of NH3. To do this, we extracted the NH3 (1,1) and (2, 2) 
spectra on a grid of points separated by 1" in the NH3 cloud 
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Fig. 10. a) 'Average' Rotational temperature map from NH3 
(1,1) and NH3 (2,2) toward MMl (see text). Scale units are 
in K. b) NH3 column density map. Scale units are in cm"^. In 
both panels red crosses mark the position of the two millimeter 
sources, VLA 8A and VLA SB, and the red tilted cross marks 
the position of IRS 1 . Note that VLA 7 lies outside the limits of 
this plot. 



MMl. We fitted the hyperfine structure of the NH3 (1,1) and 
a single Gaussian to the NH3 (2, 2) for each spectrum. For the 
NH3 (1,1) transition we fitted only the spectra with an intensity 
greater than 5 cr in order to ensure the detection of the satellite 
lines, whereas for the NH3 (2, 2) we fitted the spectra with an 
intensity greater than 4 cr. 

From the results of the fits of NH3(1,1) and NH3(2,2) 
we computed the rotational temperature (Trot) an d NH3 column 
density maps following the standard procedures (Ho & Towned 
1983; Harjuetal. 1993; Sepiilveda 1993; Angladaetal. 1991 
see Appendix for a complete description of their derivation). 
This analysis assumes implicitly that the physical conditions of 
the gas are homogeneous along the line-of-sight, i.e., the excita- 
tion and the rotational temperature are constant along the line- 
of sight. Since gradients are probably present along the line- 
of-sight (see Sect. 4.2), the values obtained from this analysis 
should be considered as some kind of 'average' along the line of 
sight. 

The map of the 'average' rotational temperature obtained is 
shown in Fig. [TOh . Interestingly, to the north of the millimeter 
sources VLA 8A and VLA 8B, there is a temperature enhance- 
ment, reaching a maximum value of 20 ± 2 K. Toward the mil- 
limeter sources the 'average' rotational temperature is around 
16 K, and it decreases toward the south. In addition, at the west- 
ern edge of the NH3 cloud MMl and toward the infrared source 
IRS 1 there is a temperature enhancement. The 'average' rota- 
tional temperature obtained at the position of IRS 1 is ~ 15± 1 K, 
which is consistent with the expected association of this source 
with the high density gas. We also found a local maximum of 
temperature around 8" to the southwest of VLA 8A. 
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Fig. 11. Top: Circularly averaged radial intensity profile at 
1 .2 mm of MM 1 . The error bars indicate the data rms inside each 
3" wide ring. The dotted lin e shows the b eam profile, including 
the error beams as given bv iGreve et al.l ()l998). The solid line 
shows the best fit model. Bottom: Spectral Energy Distribution 
(SED) for MMl in the IRAS 00213+6530 region. Filled circles 
are data from the VLA, the open circle represents the 1.2 mm 
(250 GHz) flux from the IRAM 30 m telescope, diamonds are 
IRAS data, filled diamonds are MSX data, and triangles are 
2MASS data (see Table [TTI ). The solid line is the sum of the 
free-free emission and envelope dust emission. 



FigurefTOb shows the resulting column density map for NH3, 
obtained after correction for the primary beam response. The 
highest values of the NH3 column density, 2.5x10'^ cm~^, are 
found to the south of the millimeter sources, where the rotational 
temperature shows the smallest values. Toward the two millime- 
ter sources, VLA 8A and VLA 8B, the NH3 column density is 
~ 8 X 10'"* cm-2. 

4.2. Radial intensity profile and Spectral Energy Distribution 

In order to study the spatial structure of the dusty condensa- 
tion MMl detected at 1.2 mm, we computed the circularly aver- 
aged radial intensity profile, in rings of 3" width, as a function 
of the projected distance from the 1.2 mm peak (see Table |2]). 
The result is plotted in Fig. [TT] (top panel), together with the 
IRAM 30 m beam profile, and the best model fitted to the data. 

In Table [TOl we show the parameters of the envelope mode l, 
which have been calculated by following'Estal ella et al.l (l2009h . 
In this paper the authors fit the radial intensity profile, using the 
full Planck function to describe the intensity, and the observed 
Spectral Energy Distribution (SED) simultaneously, adopting as 



Table 10. Parameters of the envelope model used to fit the radial 
intensity profile of the 1 .2 mm continuum emission and the SED 
of MMl 



parameter 


value 


Dust emissivity index yS 


1.5 


T power-law index 


0.36'' 


T at 1000 AU (K) 


31 


Density power-law index 


1.9 


Density at 1000 AU (gcm"^^) 


l.lxlO^" 


^e„v (AU) 


2.1x10* 


Envelope mass"^ (Mq) 


6.0 



" 2 1 {A +13) 

Radius for which the envelope density falls to a particle ambient 
density ^(Hz) of 7.3 x 10^ cm^^ 
Integrated up to i^env 



a model of the source a spherically symmetric envelope of gas 
and dust suiTounding the protostar(s). We do not include the 
possible contribution from the circumstellar disk because our 
angular resolution, ~ 10" (8500 AU), is much larger than the 
typical sizes of accretion disks (tens to hundreds of AU). Note 
that we do not attempt to fit the near/mid-infrared emission of 
MMl, since it originates from components at higher tempera- 
tures than that responsible for the mm/submm emission. So, the 
model is used to fit the SED up to frequencies corresponding to 
60jum. W e assumed the dust opacity la w - 0.01(v/230GHzy 
cm^ g"' (Ossenkopf & Hen nindl 19941) . being yS a free parameter 
of the model. For the density and temperature we considered 
power-laws as a function of radius, p oc r^P and T oc r~^, with p 
as a fr ee parameter of the model, and q - 2/(4+/3) (' Kenvon et al.l 
Il993h . In order to compare the model with the observed inten- 
sity profile we computed the 2-dimensional intensity map from 
the model, and we convolved it with the IRAM 30 m beam. We 
note that the beam was adopted to be the sum of two circular 
Gaussian. Then, from the convolved map we computed the cir- 
cularly averaged profile. Thus, our free parameters were the dust 
emissivity index /5, the density power-law index p, and the scale 
of the density and temperature power-laws, namely, the density 
and temperature at a radius of 1000 AU (taken arbitrarily as the 
reference radius for the power-laws). From the fitted parameters 
we derived the temperature power-law index q, the size, Rgnv, 
and mass, M^w, of the envelope. /?env is defined as the radius for 
which the envelope density falls to a particle density similar to 
the ambient density, taken as 7.3 x 10^ cm"^. Mgnv is the integral 
of the envelope mass density up to the envelope radius /?env (see 
Table[TOl). 

As mentioned above, the model fits the radial intensity pro- 
file and the SED simultaneously. The SED for MMl was built by 
using the data shown in this work and by searching the literature 
for 2MASS, IRAC-Spitzer, MSX, and IRAS data. In Table [U 
we list the photometry used for MMl (adopting the values of 
VLA 8 A for 1.3 cm, 7 mm and 2MASS, where the high angu- 
lar resolution allows us to disentangle the different sources). We 
were able to fit simultaneously the SED and the intensity pro- 
file at 1.2 mm. In Fig. [TT] (bottom panel) we show the best fit 
of the model to the observed SED. In this figure we show the 
sum of the free-free emission and the envelope flux density inte- 
grated inside the radius of the envelope. The centimeter contin- 
uum emission is dominated by free-free emission with a spectral 
index of 0.5 ± 0.4. At millimeter and submillimeter wavelengths 
the dust emission of the envelope is dominant. 
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Table 11. Photometry of the MMl clump 



1 

A 


c 

O V 


r WrlJVl 


Survey 


(jim) 


(mjy) 


( ) 


Instrument 


1.25 


< 0.2P 


2 


2MASS 


1.65 


0.431±0.00P 


2 


2MASS 


2.17 


3.104±0.001' 


2 


2MASS 


8.28 


123+10'' 


18.3 


MSX 


12 


< 250 


50 


IRAS 


12.13 


< 549 


18.3 


MSX 


14.65 


< 420 


18.3 


MSX 


21.34 


< 1200 


18.3 


MSX 


25 


980+90 


50 


IRAS 


60 


< 4910 


100 


IRAS 


100 


15000+1500 


230 


IRAS 


1200 


640+70 


11 


IRAM30m 


7000 


1.6+0.5"^ 


3 


VLA 


13000 


<0.6'='' 


1.5 


VLA 


36000 


0.27+0.02= 


12 


VLA 


60000 


0.13+0.02= 


15 


VLA 


210000 


< 2.0" 


45 


VLA' 



" J00241 110+6547095. 

G120.1368+03.0617. 

VLA8A. 
" 4(T upper limit. 

= VLA 8 (VLA 8A and VLA 8B). 

f From the NRAO VLA Sky Survey dCondon et allflgg^ . 



Table 12. Physical parameters of IRAS 00213+6530 





Mass" 


Av" 




M" 

vir 


Region 


(Mq) 


(kms-') 


(arc sec) 


(Mq) 


MMl 


3.5 


0.85 


7 


1.3 


Southern cloud 


<0.3 


0.65 


6 


0.7 


MM2 


1.1 


0.75 


11 


1.6 



" Estimated as M = /jhjMh 2 Af(H2)AA, where N(H2) is taken from 
the H2 column density map, AA is the pixel area, and a mean molec- 
ular mass per molecule /jhi = 2.8, which corresponds to a 10 % 
helium abundance. 

Intrinsic line width of the magnetic hyperfine components. 
= Deconvolved geometrical mean size of the major and minor axes 

of the source, obtained from a 2D Gaussian fit. 
" Estimated for a density distribution noc r^^. 



The model resulting from the simultaneous fit to the radial 
intensity profile and the SED is able to fit remarkably well the 
radial intensity profile, but underestimates slightly the intensity 
at projected distances ofb- 4'.'5 and 7'.'5 (see Fig[TT]top panel). 
This can be due to the presence of several YSOs close (< 10") to 
the 1.2 mm peak. The value obtained for the power-law density 
distribution index p - 1 .9 is similar to that found for other pro- 
tostellar envelopes. Chandler & Richer (2000) carried out a sub- 
millimeter survey of Class and Class I sources and fitted the 
observed radial intensity profiles with densit y index p between 
1.5 an d 2 for the majority of the sources. iHogerheiide et alj 
('2000') find p = 0.9-2 .1 in a sample of four Class I YSOs, and 
[Motte & Andrl dlOOl find p = 1.2-2.6 in a sample of embed- 
ded YSOs in Taurus and Perseus, similar to the value obtained 
here. 



4.3. Column density maps and mass 

In order to estimate the relative NH3 abundance in 100213, we 
computed the H2 column density map. To do so, we used the 
expression given in Motte et al. ( 1998) to compute the column 
density of H2 from the 1.2 mm dust in a grid of 3'.'3 x 3'.'3 (the 
pixel size of our maps), and using a dust mass opacity coeffi- 
cient K\,2mm - 0.01 cm^ g"'. Sincc the beam size of the 1.2 mm 
emission is a factor 2 larger than the beam size of the NH3 ob- 
servations, we convolved the NH3 emission to the same angular 
resolution as the 1.2 mm dust emission (~ 10"). We obtained 
the rotational temperature and column density of NH3 following 
the procedure described in Sect. 4. 1 (see Appendix). In MMl we 
obtained the rotational temperature from the NH3 (1,1) and (2,2) 
emission, which was converted to k inetic temperature using the 
relation given in lTafalla et al.l (|2004|) (see also Appendix). In or- 
der to estimate lower limits of the NH3 column density in the 
southern cloud and MM2 we adopted a rotational temperature 
of 10 K (see Sect. 5.2). 

Maps of the H2 column density and NH3 column density are 
shown in Fig. [12] The H2 column density has a maximum value 
of ~ 4 X 10^^ cm"^ toward the peak position of the 1 .2 mm emis- 
sion of MMl, and decreases in the more extended structure to 
values of ~ 10^^ cm"^. The uncertainty in the H2 column den- 
sity is around a factor of 2, mainly due to the uncertainty in the 
dust opacity, as the uncertainty in the mm flux density (the rms 
of the map), and in the dust temperature constitute only a small 
contribution to the total uncertainty. The NH3 column density 
map, corrected for the primary beam response, also shows differ- 
ences between the three NH3 clouds (see Fig[T2b). In MMl, the 
maximum value of the NH3 column density, ~ 1.5 x lO''' cm"^, 
is found to the south of VLA 8A and VLA SB, which is consis- 
tent with the results found in Sect. 4.1. The uncertainty in the 
NH3 column density is of the order of 10 - 20 %. It has been 
estimated from the uncertainty in rotational temperature, in line 
width, and in At,,, (see Appendix). The values of the NH3 col- 
umn density toward MM2 are significantly higher, in the range 
of (2-7)xl0'''. For the case of MM2 and the southern cloud, the 
NH3 column density derived is a lower limit since the rotational 
temperature adopted is an upper limit due to the non-detection 
ofNH3 (2,2). 

We estimated the mass of each condensation from the H2 
column density maps (see Table \T2\ . The mass in MMl is 
~3.5 Mq, consistent with the values obtained from the 7 mm 
continuum emission and for the envelope model (~ 6 Mq). In 
addition, we compared the valu es obtained wi t h the virial mass, 
Mvir, estimated using Eq. (5) of iBeltran et al.l (l2006h . which as- 
sumes a spherical cloud with a power-law density distribution 
p oc r^P, with p= 2.0, and neglecting contributions from mag- 
netic fields and surface pressure. As can be seen in Table [12] 
the mass of MMl is higher than the virial mass indicating that it 
might be unstable and undergoing collapse. The mass of MM2 
is M ^ 0.7 Mvir, indicating that the material in this cloud is 
stable. In contrast, toward the southern cloud the total mass 
of gas is M < 0.4 Mvir, being this clump gravitationally un- 
bound and could disperse at roughly the internal sound speed of 
~ 0.3 km s"' on a timescale of around 10^ yr, unless it is con- 
fined by external pressure. All this suggests that it could be in an 
earlier evolutionary stage than MMl and MM2. 

4.4. Relative NH3 abundance 

In Fig. [T2t we show the relative NH3 abundance maps of MMl 
and MM2. The typical value of the NH3 abundance found in 
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Fig. 12. a) H2 column density map from the 1.2 mm dust emission, b) NH3 column density map. c) NHs/Hiabundance map, for 
MMl and the southern cloud (left panels) and for MM2 (right panels) in logarithmic scale. Scale units for the H2 and NH3 column 
density are cm"^. In all panels red crosses indicate the position of the two millimeter sources, VLA8A and VLA8B, and the red 
tilted cross marks the position of the infrared source IRS 1. The blue cross denotes the position of VLA 7. The color scale is the 
same for the left and right panels. 



MMl, in a position just south of the three YSOs, is around 
2 X 10"^ - 4 X 1 0"^, which is similar to the typical value found 
in dense clouds jHerbst & Klempere3ll973L see Anglada et al. 
1995 for a discussion on NH3 abundances). In MMl, near the 
YSOs, there is a slight trend in the NH3 abundance to increase, 
from 1 X 10 near IRS 1, up to 3 x 10 toward VLA 8B. In 
this case, the most evolved YSO IRS 1 is likely dispersing the 
dense gas material of its suiToundings, and the result is a de- 
crease of the NH3 abundance. It is, however, a modest effect, 
slightly above the typical uncertainty, of a factor of 2. Toward 
the southern cloud we derived a lower limit for the relative NH3 
abundance of ~ (1 - 4) x 10"^. 

Regarding the NH3 abundance in MM2, we find values rang- 
ing from 1 - 2 X 10"^ up to 1 x 10"^. This abundance is higher 
than in MMl, despite the observed NH3 emission in MM2 be- 
ing fainter than in MMl, mainly due to two effects: the NH3 
column density is higher due to the correction for the primary 
beam response, and the H2 column density is low, ~ 3 x 10^^ - 
1 X 10^^ cm"^. It is important to remark that the relative NH3 



abundan ce in MM2 is a lower lim it. Al though these abundances 
are high. lBenson & MversI (Il983h and lOhishi elal] (Il992h find 
NH3 abundances around ~ 2 x 10"^ in starless cores of the dark 
clouds L1498 and L1512, and in L134N,_respectively. In addi- 
tion, chemical models of iHartquist. Williams. & Vitil ( I2b0ll) are 
able to reproduce such a high NH3 abundance for a young core. 
Core D in TMC-1, with N2 and CO not freezing-out. Therefore, 
it is plausible that MM2 and the southern cloud are in fact evo- 
lutionary young, as suggested by the starless properties of this 
cloud (see Sect. 5.2). 



5. Discussion 

The results obtained for IRAS 00213+6530 show that in this 
region there are different sources, which have different radio and 
infrared properties: 
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5.1. VLA 8A, VLA 8B, and IRS1: Multiple sources in different 
evolutionary stages 

The centimeter and millimeter continuum observations, together 
with near-infrared data, allowed us to identify three YSOs in 
MMl: IRS 1, VLA 8A and VLA 8B. While at 3.6 cm we could 
not resolve the centimeter emission of VLA 8, at 7 mm we found 
two sources VLA 8 A and VLA 8B. In the near infrared, only 
VLA 8A and IRS 1 are detected. In addition, all sources are 
deeply embedded in the dense gas traced by NH3. Line broad- 
ening and local heating have been detected toward their posi- 
tion, indicating the true association of the three objects with the 
molecular gas. 

From these preliminary results we can make a rough estimate 
of the evolutionary stage of the detected sources. First, IRS 1 is 
associated with an infrared source with no detected 7 mm con- 
tinuum emission. VLA 8 A, although bright in the infrared, is the 
strongest millimeter source in the field. Its SED has a quite steep 
profile at the 2MASS wavelengths. Thus, d(log/lF^)/d(log/l) > 
between 1 and 10 yum, which is consistent with the classification 
of VLA 8A as a Class I source (e. g., Hartmann 1998). Regarding 
the near-infrared colors of the 2MASS sources associated with 
IRS 1 and VLA 8A, derived from 2MASS photometry (see 
Table [T3Tl. we found that both sourc es fall inside the area corre- 
sponding to YSOs of Class O/I (e. g..' Itoh et al.lll996tlOiha et al.1 
[2004; Matsuyanagi et al. 2006). Finally, VLA 8B, which is also 
associated with dense gas tracers, shows no infrared emission 
at all. At 7 mm VLA 8B has little dust continuum emission 
associated, suggesting that the dust emission is colder, proba- 
bly detectable at 1 mm and/or submillimeter wavelengths. Thus, 
VLA 8B could be in a previous stage of evolution, being still 
more embedded than VLA 8A, possibly in the Class phase. It 
is worth noting that we found an extended temperature enhance- 
ment to the north of VLA 8A and VLA 8B, which seems to be 
associated with the passage of an outflow that heats and perturbs 
the dense gas: Bus quet et al. (in pre p.) find that the large scale 
molecular outflow (I Yang et al.ll 19901) . when observed with high 
angular resolution, is centered on a position near VLA 8A and 
VLA 8B, being these sources the candidates to be driving the 
observed high-velocity gas. 

Thus, it seems that the 100213 region is harboring a multi- 
ple system of low mass protostars, indicating that the star for- 
mation process in this region is not producing a single YSO. 
Since low mass protostars evolve approximately at the same 
rate to the main sequence, the different evolutionary stages 
found in 100213 suggest that stars in this region are not form- 
ing simultaneously but continuously. Actually, there may be 
different generations due to different timescales of core col- 
lapse, as has been f ound in other low mass star-fo rming re- 
gions (e.g., L1551: iMoriartv-Schieven et al.l l2006h . indicat- 
ing that the formation of different stars is not simultaneous 
but sequential in time, maybe triggered by the interaction of 
the molecular outflow with a dense core in its surroundings 
jYokogawa et al.ll2003l:IShimaiiri et al ]l2008l) . Therefore, the ini- 
tial assumption that star formation occurs in an isolated mode 
may not be appropriate to describe the 100213 region when 
the region is studied with high angular resolution. This poses 
the question of what extent can we adopt the isolated mode 
in the theories of star formation, as this high angular res- 
olution s tudy together wit h a large number of recent stud- 
ies (e. g., Huard et al. 1999; Diupvik et al. '2006; Teixeira et al.' 
2007; Carrasco-Gonzalez et al. 2008; Chen et al. 2008, 2009; 
Forbrich et al. 2009; Gut ermuth et all l2008t ISwift & WelchI 
200&: iGirart et al.. ,2009.) suggest that isolated star formation 



Table 13. Infrared excess of 2MASS sources 



2MASS 


Associated 








Source 


Source 






(H-K) 


J00241110+6547095 


VLA 8A 


1.02d 


:0.16 


1.47±0.18 


J00241010+6547091 


IRS 1 


1.25d 


:0.26 


2.61±0.27 


J00241251+6546418 




0.48d 


:0.04 


0.12±0.05 



seems to be rare in the Galaxy, even in low mass star-forming 
regions. 

It is interesting to note that the three YSOs are spatially or- 
dered from youngest (east) to oldest (west), suggesting that an 
external agent could be inducing star formation in MMl. For 
this, it would be very useful to identify and map the molecular 
outflows in the region. Alternatively, it would be very useful as 
well to study the possible association of 100213 with the HII re- 
gion S 1 7 1 , located to the northwest of 1002 1 3 . Finally, it is worth 
noting that the 1002 1 3 region falls exactly on the southern border 
of the Cep OB4 shell (lKunll2008l) . 

5.2. Starless candidates: MM2 and the southern cloud 

While the southern cloud is detected only in NH3, MM2 is 
detected both in NH3 and dust emission. Given that the (2, 2) 
line is not detected toward southern cloud, and only detected 
toward the northern peak of MM2, these clouds are cold, 
since r,ot <11 K. As mentioned in Sect. 3.3, the near-infrared 
source 2MASS J0024125 1-1-6546418 spatially coinciding with 
the southern cloud is not likely associated with the dense gas. We 
estimated the infrared excess from the ( J-H) vs (H-J) diagram 
(see Table [T3]). The near-infrared colors derived are characteris- 
tic of main sequence stars, giants, supergiants. Class III sources , 
or Class II sources with small infrared excess (Itoh etalJ [T995 
iMatsuvanagi et al]|2006l) . In addition, this source has an optical 
counterpart seen in the DSS2 image. Thus, it is likely a fore- 
ground source, not associated with the NH3 dense gas and the 
100213 star-forming region. Regarding MM2, no near-infrared 
sources are associated with it. In addition, neither MM2 nor the 
southern cloud seem to be associated with molecular outflows 
(Busquet et al. in prep.). Therefore, no clear signposts of stellar 
activity are found for these two clouds, suggesting they could be 
starless. 

5.3. On the nature of VLA 7 

VLA 7, lying outside the NH3 dense gas and dust emission, 
has a very negative spectral index at centimeter wavelengths 
(-1.6 ± 0.2), indicating that the emission has a non-thermal ori- 
gin, found typically for extragalacti c sources with steep spec- 
trum and some pulsars (e. g.. ILehtinen et al.ll2003h . Given the 
close proximity to MMl, we considered the possibility that 
VLA 7 could be related to the 100213 region. In low mass star- 
forming regions non-thermal emission has been detected toward 
some YSOs , like Class YSOs (Choietal. 2008) or T-Tauri 
stars (e. g., lAndrel 119961; iRodrfguez et al.1 1 19991: iGibbI Il999h . 
However, centimeter emission arising from weak T-Tauri stars 
is usually polarized, and they are often optically visible. We did 
not find evidences of circular polarized emission toward VLA 7, 
neither a visible counterpart, so we consider unlikely this pos- 
sibility. Another possibility to explain the negative spectral in- 
dex of VLA 7 is non-thermal synchrotron emission produced 
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in shocked regions of outflowin g gas, found mainly iri high 
mass star-forming regions (e. g.. Rodriguez & Reipurthl ll989t: 
[Oaray et al. 2003; Rodriguez et al. 2005). As there is molecu- 
lar outflow emission in this region (Busquet et al. in prep.) we 
can not rule out this possibility, and further observations would 
help to confirm the truly association of VLA 7 with the 100213 
star-forming region. 

6. Conclusions 

We observed with the VLA, IRAM30 m Telescope, and the 
NASA 70 m antenna the continuum emission at 6 cm, 3.6 cm, 
1.3 cm, 7 mm, and 1.2 mm, the NH3 (1,1) and NH3 (2,2) lines, 
and the H2O maser and CCS emission toward the low mass star- 
forming region IRAS 00213+6530. Our main conclusions can 
be summarized as follows: 

1. The 1.2 mm continuum emission observed with the 
IRAM 30 m shows two dust condensations, MMl and MM2. 
The continuum emission at centimeter and millimeter wave- 
lengths, together with the available data from 2MASS, have 
revealed three sources, IRS 1, VLA 8 A, and VLA SB, all 
embedded in the dusty cloud MMl. These sources show dif- 
ferent radio and infrared properties, and seem to be in differ- 
ent evolutionary stages, with VLA 8B being in the earliest 
phase. In MMl, low mass star formation appears to proceed 
along a west-east direction. 

2. We marginally detected H2O maser emission toward 100213 
with the NASA 70 m antenna during the observations caiTied 
out on 2008 September 23, but other attempts, including the 
VLA observations, yielded negative results. 

3. The YSOs found in the region are deeply embedded in the 
high-density gas. The NH3(1,1) emission traces an elon- 
gated structure that consists of a main cloud (MMl and the 
southern cloud) and a smaller cloud, MM2, located to the 
northwest of MMl. While the southern cloud and MM2 ap- 
pear as quiescent and starless, in MMl there is evidence of 
a perturbation of the gas (line broadening and local heat- 
ing) along the east-west direction, associated with IRS 1 and 
north of the two millimeter sources, elongated in the north- 
south direction. We propose that part of the dense gas is be- 
ing perturbed by the passage of one or more outflow(s). 

4. We detected CCS emission toward 100213 using the 
NASA 70 m antenna with a line width, ~ 1.3 km s"', large 
compared with previous studies. 

5. The source VLA 7, which has a negative spectral index, lies 
outside but near the border of the NH3 (1,1) condensation. 
Although this source could be a background source we can 
not rule out the possibility that VLA 7 could be in fact the 
result of the interaction of a molecular outflow with the sur- 
rounding medium. 

6. We used a spherically symmetric envelope model that fits 
simultaneously the observed SED from 7 mm to 60 jum, and 
the radial intensity profile at 1.2 mm of the clump associated 
with MMl. The best fit was obtained for a dust opacity law 
index fi - 1.5, a temperature at 1000 AU of 31 K, and a 
density at 1 000 AU ofl.lxlO"''^g cm"-' or a particle density 
of «(H2)~ 2.3 X lO'* cm""*. The envelope radius is /?env ~ 
21000 AU (3.1x10'^ cm), and inside this radius the envelope 
model mass is 6 Mq. 

7. There is a strong differentiation of NH3 abundance in the 
region. In particular, we found low values, ~ 2 x 10"^ of 
the NH3 abundance associated with MMl, which contains 
the YSOs. On the other hand, toward those clouds with star- 



less properties (the southern cloud and MM2) the NH3 abun- 
dance rises up to ~ 1 x 10"^, suggesting that in evolved clouds 
with star-formation activity there is a decrease in the NH3 
abundance. 

Acknowledgements. G.B. is grateful to Serena Viti, David A. Williams, and 
Oscar Morata for useful discussion on the ammonia abundance. We are grate- 
ful to an anonymous referee and to the editor for valuable comments. The au- 
thors are supported by the Spanish MEC grant AYA2005-08523-C03, and the 
MICINN grant AYA2008-06189-C03 (co-funded with FEDER funds). A.P. is 
also supported by the MICINN grant ESP2007-65475-C02-02 and the program 
ASTRID S0505/ESP-0361 from La Comunidad de Madrid and the European 
Social Fund. G.A. acknowledges support from Junta de Andalucfa. This pub- 
lication makes use of the data products from the Two Micron All Sky Survey, 
which is a joint project of the University of Massachusetts and the Infrared 
Processing and Analysis Center/CaHfornia Institute of Technology, funded by 
the National Aeronautics and Space Administration (NASA) and the National 
Science Foundation. 



References 

Adams, F C, & Myers, P. C. 2001, ApJ, 553, 744 

Andre, R 1996, in Taylor A.R., Paredes, J.M., eds, ASP Conf Ser Vol 93, Radio 

Emission from the Stars and the Sun. Astrom. Soc. Pac, San Francisco, p. 273 
Anglada, G., Estalella, R., Mauersberger, J., Torrelles, J. M., Rodriguez, L. F, 

Canto, J., Ho, R T. R, & D'Alessio, R 1995, ApJ, 443, 682 
Anglada, G., Sepulveda, I., & Gomez, J. F 1997, A&AS, 121, 255 
Anglada, G., Villuendas, E., Estalella, R., Beltran, M. T, Rodriguez, L. F, 

Torrelles, J. M. & Curiel, S. 1998, ApJ, 116, 2953 
Beltran, M. T., Estalella, R., Anglada, G., Rodriguez, L. E & Torrelles, J. M. 

2001, ApJ, 121, 1556 
Beltran, M. T, Girart, J. M., & Estalella, R. 2006, A&A, 457, 865 
Benson, P J., & Myers, P C. 1983 ApJ, 270, 589 
Benson, P J., & Myers, P C. 1989 ApJS, 71, 89 

Briggs, D. 1995, PhD Thesis, New Mexico Inst, of Mining and Technology 
Brooke, T. Y., Huard, T. L., Bourke, T. L., Boogert, A. C. A., Allen, L. E., Blake, 

G. A., Evans, N. J., II, Harvey, P. ,M., Koerner, D. W., Mundy, L. G. et al. 

2007, ApJ, 655, 364 
Carilli, C. L., & Holdaway, M. A. 1997, Millimeter Array Technical Memo 173 

(NRAO) 

CaiTasco-Gonzalez, C, Anglada, G., Rodriguez, L. E, Torrelles, J. M., Osoiio, 

M., & Girart, J. M. 2008, ApJ, 676, 1073 
Chandler, C. J., & Richer, J. S. 2000, ApJ, 530, 851 

Chen, X., Bourke, T. L., Launhardt, R., & Henning, T 2008, ApJ, 686, L107 
Chen, X., Launhardt, R., & Henning, T. 2009, ApJ, 691, 1729 
Choi, M., Hamaguchi, K., Lee, J-.E., & Tatematsu, K. 2008, ApJ, 687, 406 
Clarke, C. J., Bonnell, I. A., & Hillenbrand. L. A. 2000, in Protostars £r Planets 

IV, eds. V. Mannings, A. P. Boss, & S. S. Russell, (Tucson: University of 

Arizona Press), p. 151 
Condon, J. J., Cotton, W. D., Greisen, E. W., Yin, Q. E, Periey, R. A., Taylor, 

G. B., & Broderick, J. J. 1998 AJ, 115, 1693 
de Gregorio-Monsalvo, I., Chandler, C. J., Gomez, J. F, Kuiper, T. B. H., 

Torrelles, J. M., Anglada, G. 2005, ApJ, 628, 789 
de Gregorio-Monsalvo, I., Gomez, J. E., Suarez, O., Kuiper, T. B. H., Rodriguez, 

L. F, Jimenez-Bailon, E. 2006, ApJ, 642, 319 
Djupvik, A. A., Andre, Ph., Bontemps, S., Motte, E, Olofsson, G., Galfalk, M., 

Floren, H.-G. 2006, A&A, 458, 789 
Estalella, R., Patau, A., Girart, J. M., Beltran, M. T, Osorio, M., Ho, P T. P & 

Anglada, G. 2009, A&A, submitted 
Eelli, M., Palagi, E, & Tofani, G. 1992, A&A, 255, 293 

Eorbrich, J., Stanke, T, Klein, R., Henning, T., Menten, K. M., Schreyer, K., 

Rosselt, B. 2009, A&A, 493, 547 
Gai-ay G., Brooks K. J., Mardones D., & Nonis R. P 2003, ApJ, 587, 739 
Gibb, A. G. 1999, MNRAS, 304, 1 

Girart, J. M., Estalella, R., Ho, P T. P, & Rudolph, A. L. 2000, ApJ, 539, 763 

Girart, J. M., Rao, R., & Estalella, R. 2009, ApJ, 694, 56 

Gomez, M., Hartmann, L., Kenyon, S. J., & Hewett, R. 1993, AJ, 105, 1927 

Greve, A., Kramer, C, & Wild, W. 1998, A&AS, 133, 271 

Gutermuth, R. A., Bourke, T. L., Allen, L. E., Myers, R C, Megeath, S. T, 

Matthews, B. C. et al. 2008, ApJ, 673, L151 
Han, E, Mao, R. Q., Lu, J., Wu, Y. E, Sun, J., Wang, J. S., Pel, C. C, Fan, Y, 

Tang, G. S. & Ji, H. R. 1998, A&AS, 127 181 
Harju, J., Walmsley, C. M., & Wouterioot, J. G. A. 1993, A&AS, 98, 51 
Hartmann, L. 1998, in Accretion Processes in Star Formation, eds. A. King, 

D. Lin, S. Maran, J. Pringle, & M. Ward, Cambrigde: Cambridge University 

Press, p. 8 



16 



G. Busquet et al.: Star formation in IRAS 00213+6530 



Hartquist, T. W., Williams, D. A., & Viti, S. 2001, A&A, 369, 605 

Herbst, E., Klemperer, W. 1973, ApJ, 185, 505 

Ho, P. T. P., & Townes, C. H. 1983, ARA&A, 21, 239 

Hogerheijde, M. R. & Sandell, G. 2000, ApJ, 534, 880 

Hotzel, S., Harju, J., & Walmsley, C. M. 2004, A&A, 415, 1065 

Huard, T. L., Sandell, G., & Weintraub, D. A. 1999, ApJ, 526, 833 

Itoh, Y., Tamura, M., & Gatley, I. 1996, ApJ, 465, L129 

Kenyon, S. J., Calvet, N., & Hartmann, L. 1993, ApJ, 414, 676 

Kun, M. 2008, Star Forming Regions in Cassiopeia Handbook of Star Forming 

Regions, Volume I: The Northern Sky ASP Monograph Publications, 

Vol. 4. Edited by Bo Reipurth, p.240 
Lada, C. J. 1999, in The Origin of Stars and Planetary Systems, eds. C. J. Lada, 

and N. D. Kylafis, Kluwer Acad. Publ., p. 143 
Lada, C. J., & Lada, E. A. 2003, ARA&A, 41, 57 

Lada, E. A., Strom, K. M., & Myers, P C. 1993, in Protostars &■ Planets III, eds. 

E. H. Levy and J. 1. Lumine, Tucson: University of Arizona Press p. 245 
Lee, J-E., Di Francesco, J., Lai, S-P, Bourke, T. L., Evans 11, N. J., Spiesman, B., 

Myers, P C, Allen, L. E., Brooke, T. Y, Porras, A., & Wahhaj, Z. 2006, ApJ, 

648, 49 

Lehtinen, K., Haiju, J., Kontinen, S., & Higdon, J. L. 2003, A&A, 401, 1017 
Matsuyanagi, I., Itoh, Y, Sugitani. K., et al. 2006, PASJ, 58, L29 
Moriarty-Schieven, G. H., lohnstone, D., Bally, J., & Jenness, T. 2006, ApJ, 
645,357 

Motte, F, & Andre, P 2001, A&A, 365, 440 
Motte, F, Andre, P, & Neri, R. 1998, A&A, 336, 150 

Ohishi, M., Irvine, W. M., & Kaifu, N. 1992, in Astrochemistry of Cosmic 

Phenomena, ed. PD. Singh (Dordrecht: Kluwer), 171 
Ojha, D. K., Tamura, M., Nakajima, Y, & Fukagawa M. 2004, ApJ, 608, 797 
Ossenkopf V, & Henning, Th. 1994, A&A, 291, 943 

Pauls, A., Wilson, T. L., Bieging, J. H., & Martin, R. N. 1983, A&A, 124, 23 
Pfalzner, S., Tackenberg, J., & Steinhausen, M. 2008, A&A, 487, L45 
Rodriguez, L. F, Anglada, G. & Curiel, S. 1999, ApJS, 125, 427 
Rodriguez L. F, Garay G., Brooks K. J.,& Mardones D. 2005, ApJ, 626, 953 
Rodriguez L. F, & Reipurth, B. 1989, Rev. Mex. Astron. Astrofis. 17, 59 
Sepulveda, 1. 1993, Master Thesis. University of Barcelona 
Sepiilveda, 1. 2001, PhD Thesis. Universitat de Barcelona 
Shimajiri, Y, Takahashi, S., Takakuwa, S., Saito, M., & Kawabe, R. 2008, ApJ, 
683, 255 

Shu, F H., Adams, F C, & Lizano, S. 1987, ARA&A, 25, 23 

Sweitzer J. S. 1978, ApJ, 225, 116 

Swift, J. J., & Welch, W. J. 2008, ApJS, 174, 202 

Tafalla, M., Myers, P C, Caselli, P, & Walmsley, C. M. 2004, A&A, 416, 191 
Teixeira, P S., Zapata, L. A., & Lada, C. J. 2007, Apl, 667, L179 
Ungerechts, H., Winnewisser, G., & Walmsley, C. M. 1986, A&A, 157, 207 
Yang, J., Fukui, Y, Umemoto, T. & Ogawa, H. 1990, ApJ, 362, 538 
Yokogawa, S., Kitamura, Y, Momose, M., & Kawabe, R. 2003, ApJ, 595,266 



Appendix A: Derivation of Trot and A'^(NH3) from 
NH3(1,1) and (2,2) observations 

rMB(l,l;'w) and A^(l,l); 

The NH3(1,1) method of CLASS fits the magnetic hyperfine 
structure of NH3 (1,1). The output parameters for the fit to the 
hyperfine structure are: At^,, the velocity of the reference line, 
the intrinsic line width, and the optical depth of the (1,1) main 
line, Tm (sum of the optical depths of the magnetic hyperfine 
components of the main lin e, - t(1, 1)/2 . The parameter A 
is defined as, according to Pauls et al] (Il983h . A = /[A(7'ex) - 
JyiJb^}, where / is the filling factor Then, from the output pa- 
rameters, and applying the radiative transfer equation, one can 
obtain the (1,1) main line temperature Tmb (1,1; m). 



1 



(A.l) 



The excitation temperature Tex is not directly obtained from 
the fit, but is isolated from the output parameter Ar^, 



T - 

^ ex — 



1.14 



ln(l + 1.14/[AWT„ + A,(rbg)])' 



Note that no assumption is made concerning Tex with respect to 
the background temperature Tbg. 

The beam aver aged column density in the (1,1) level 
dAnglada etalj|1995h . 



A^(l,l) 



-,-2 



= 1.58X 10'7-r 



14/7-„ _ I 



Av 



km s 



(A.3) 



the filling factor / being assumed to be 1 for our VLA observa- 
tions. 

To derive Eq. (IA.3b . A^( 1 , 1 ) is not approximated to 2N+ (1,1), 
but is taken as A^(l, 1) = A^+(l, 1)[1 -i-exp(/jvii/^rex)] (see Harju 
et al. (1993) for more details). 

Tub (2,2;m): 

For NH3 (2,2) we fitted one single Gaussian, being the (2,2) main 
line temperature, Tmb (2, 2; m), an output parameter of the fit. 



n21 . 



The rotational temperature derived from NH3(1,1) and 
NH3 (2,2) can be estimated, following Ho & Townes (1983, 
Eq. 4), by assuming that the transitions between the metastable 
inversion doublets are approximated as a two-level system, and 
that the excitation temperature Tex and line width Av are the 
same for both NH3 (1,1) and NH3 (2,2). Then, 



n21 



-41.5 



In 



(_o^ln[l 



Tub (2,2;m) 
Tub (l,l;ni) 



(1 



-.1) 



(A.4) 



Note that we did not assume that the emission is optically thin. 
The assumption of a two-level system is reasonable because 
transitions between the metastable inversion doub lets are usually 
much faster than those to other rotational states (iHo & To wnesI 
Il983h . If the density and temperature were high enough to popu- 
late the upper non metastable st ates, multilevel statistical calcu- 
lations would be required (e. g.. lSweitzer et al.| [T978). 

An estimate of the gas kinetic temperature can be obtained 
by correcting the rotational temperature derived from NH3, using 
the expression given in Tafalla et al. (2004), 



n21 



l-glnfl + l.le-if'/T-.l',]' 



(A.5) 



which is almost independent of core density and size. This rela- 
tion is recommended for the range Tk - 5-20 K. 

N(NH3): 

The NH3 column d ensity was derived by following 
lUngerechts eTal] (11986 !). and Harju et al. (1993). The main 
assumptions are: ;) only metastable levels are populated; //) 
Ti-ot is the same for each pair of rotational levels; the ratio 
of the column densities of each rotational level is the same 
to the ratio of the column densities of upper inversion levels; 
iv) the contribution to the total NH3 column density comes 
essentially from levels with J < 3; v) the relative population of 
all metastable levels of both ortho and para-NH3 is that given by 
thermal equilibrium at temperature Trot; and vi) the frequencies 
for the NH3(1,1) and NH3 (2,2) transitions are very similar. 
With these assumptions. 



(A.2) A^(NH3) =A^(1,1) 



1 23 4/7"-' 1 ^ 
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Uncertainty in T^^^: 

In order to estimate the uncertainties associated with Trot and 
A^(NH3) introduced by this method, we did the following. 

The error of was estimated by assuming opti- 
cally thin emission and that the main sources of er- 
ror come from rMB(l,l;w) and rMB(2,2;m). Defining 
R = 7'mb (2, 2;ot)/7'mb(1, liw), the relative error is 

f = ^|{'-m^f H^^§Mf^ <^rMB(U;m) and 
STmb (2, 2; m) given directly by the hyperfine fit. Then, the er- 
ror in the rotational temperature has been estimated as, 

„ -41.5 SR 
ln^(0.283/?) R 

Then, as a test for the previous estimate of the error in Trot, 
we estimated the opacity from the ratio of the main line intensity 
to the inner satellites averaged intensity, Tmb (1,1, is), following 
Ho & Townes (1983), 

TMB(l,l;Ji) l-e-W3.6' 

and derived T^ot with this estimate of the opacity. The opacities 
derived from this method are systematically lower but compati- 
ble with the values derived from the hyperfine fit, and the rota- 
tional temperature obtained agrees with the values derived from 
the NH3 (1,1) hyperfine method used in this work. We note that 
the opacity inferred from the hyperfine fit seems to be more reli- 
able than the opacity from the ratio of the main line to the satel- 
lites when the width of the magnetic hyperfine components is 
comparable to the their separation in velocity (e. g., Anglada et 
al. 1995), which is 0. 1 1-0.53 km s ' . This is probably the case of 
our region. However, since the observations reported here were 
carried out with a spectral resolution in some cases lower than 
the intrinsic line width, the hyperfine fits must be regarded with 
caution and for this reason we compared them with the ratio of 
the main line to the satellites method. 



